Radiative heat transfer between dielectric bodies 
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The recent development of a scanning thermal microscope (SThM) has led to measurements 
of radiative heat transfer between a heated sensor and a cooled sample down to the nanometer 
range. This allows for comparision of the known theoretical description of radiative heat transfer, 
which is based on fluctuating electrodynamics, with experiment. The theory itself is a macroscopic 
theory, which can be expected to break down at distances much smaller than 10 _8 m. Against 
this background it seems to be reasonable to revisit the known macroscopic theory of fluctuating 
electrodynamics and of radiative heat transfer. 

PACS numbers: 44.40.+a, 05.40.-a, 03.50.De 



o 
o 



> 

(N 
cn 

rn 

o 



X 



The development of a new kind of scanning thermal mi- 
croscope (SThM) by the experimental Energy and Semi- 
conductor Research group in Oldenburg has made possi- 
ble measurements of radiative heat transfer down to the 
nanometer range. A heated sensor tip incorporating a 
thermo-couple with a tip radius of 60nm scans a cooled 
sample surface, registering the thermo- voltage in the sen- 
sor. In order to assure that the results are not influenced 
by convection or adsorbates on the sample surface, the 
SThM works under ultra high vacuum conditions. In 
the future, scanning thermal microscopy might become 
a technique supplementing the well-known scanning tun- 
neling and atomic force microscopy. On the other hand, 
the measurements with the SThM allow for comparision 
of the existing theory of radiative heat transfer with the 
experiment. This shall be motivation enough for revisit- 
ing the theory of radiative heat transfer. 

INTRODUCTION 

The radiative heat transfer between macroscopic 
dielectric bodies is a well-understood physical phe- 
nomenon! and the energy transported per unit time 
and unit area can be expressed by the Kirchhoff-Planck 
law [l[ in the stationary case, i.e., when the bodies are in 
local thermal equilibrium. But this is only true for dis- 
tances d which are large compared with the thermal wave 
length A t h given by Wien's law, i.e., for d 3> Ath- For 
distances d <C A t h the radiative heat transfer increases 
strongly due to so-called evanescent modes. 

For convenience we consider two semi-infinite bodies 
with a vacuum gap of distance d between them (see 
Fig. [IJ. This system is invariant under rotations and 
we chose the z-axis as the axis of rotation. We further 
assume that the first body at z < has a temperature 
Ti =/= 0, whereas the second semi-infinite body at z > d 
has zero temperature, i.e., Ti = 0. Due to quantum and 
vacuum fluctuations of the atoms and electrons in the 
first body there will be a fluctuating electromagnetical 
field inside that body. Now, the boundary conditions of 
macroscopic electrodynamics require that the transver- 




FIG. 1: Sketch of the propagating and evanescent modes con- 
tributing to the radiative heat transfer between two semi- 
infinite bodies at different temperatures. 



sal components of the electric and the magnetic field 
are continuous at the interface z = 0. Therefore, there 
will also be a fluctuating electromagnetic field outside 
the semi-infinite body produced by fluctuating internal 
sources. This field can be divided into a propagating and 
an evanescent part. This becomes clear if we consider a 
plane wave solution exp(i(k z z — Lot)) in the vacuum gap 
propagating into the positive z-direction, with the mod- 
ulus of the wave vector given by 




(1) 



Obviously, the modulus of the wave vector k z is purely 
real if the modulus of the transversal wave vector k± is 
smaller than the modulus of the vacuum wave vector, 



In this case we will have an oscillatory solu- 



tion, i.e., propagating waves. On the other hand, k z is 
purely imaginary if > fco ■ In that case the plane wave 
is exponentially damped in z-direction, i.e., the solution 
is a so-called evanescent wave. 

The characteristic wave length of the fluctuating field 
in the vacuum gap, generated in the semi-infinite body 
at z < 0, is the thermal wave length Ath- It follows 
that for distances d 3> A, i.e., in the far-field region, the 
evanescent field can not reach into the second body and 
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therefore only propagating modes can contribute to the 
radiative energy transfer. This energy transfer — as men- 
tioned above — can then be described with the Kirchhoff- 
Planck law and is smaller than the energy transfer be- 
tween two black bodies, which is given by the Stefan- 
Boltzmann law [l[ 



Sbb = &{Ti - T 2 4 ) 



(2) 



with a = 5.67 • 10" 8 Wm" 2 K" 4 . On the other hand, for 
distances d between the two bodies smaller than Ath, i-e., 
in the near-held region, the evanescent waves generated in 
the first body can reach into the second body and prop- 
agate there. This can be seen as some kind of photon 
tunneling through a vacuum barrier (see Fig. [1]). Hence, 
in the near-field the propagating and evanescent fluctu- 
ating fields both will contribute to the radiative energy 
transfer. In that case, the energy transfer can be several 
orders of magnitude larger than the black body value ((2|) , 
as we will see later. 



THEORETICAL DESCRIPTION 
Fluctuating fields 

The theoretical treatment of radiative heat transfer is 
usually based on Rytov's fluctuating electrodynamics Q ■ 
Within this framework the macroscopic Maxwell equa- 
tions are augmented by fluctuating source currents j , de- 
scribing the fluctuating sources of the electric and mag- 
netic field E and H inside the dielectric dissipative body. 
The individual frequency components of these source cur- 
rents j(r, uj) are considered as stochastic Gaussian vari- 
ables. Within this treatment the Maxwell equations be- 
come so-called stochastic Maxwell equations 

V x E(r,cj) = icj^ H(r,w), 
VxH(r,w)=j(r,w)-iw£(u)E(r,w) (3) 

with the permeability of the vacuum fio and the per- 
mittivity of the body e(w). By writing the stochastic 
Maxwell equations in the form of eq. (J3|) the dielectric 
body containing the source currents is already assumed 
to be non-magnetic, homogeneous, isotropic and local. 

The electric and magnetic fields described by the 
stochastic Maxwell equations are classical stochastic pro- 
cesses and can formally be expressed as 

E(r,u) = iujfx /dV G £ (r,r')j(r», 

H(r, W )=k^ Jd 3 r'G H (r,r')i(r',u) (4) 

where the integral has to be taken over the source region, 
i.e., over the volume of the dielectric body containing the 
source currents. The tensors G E and G H are the dyadic 
electric and the magnetic Green's function. The electric 



dyadic Green's function is a solution of the Helmholtz 
wave equation 



(V x V x -fc 2 )G £ (r,r') = 16(r - r') 



(5) 



with the wave vector k and the unit matrix 1. By def- 
inition, the electric dyadic Green's function satisfies the 
boundary conditions of the electric field. Even though a 
similiar statement holds for the magnetic Green's func- 
tion, it is more convenient to construct the dyadic mag- 
netic Green's from the electric Green's function by means 
of the relation 

G H (r, r') = — V x G £ (r, r'), (6) 

which is a direct consequence of Maxwell's equations. 

The reformulation of the electric and magnetic field E 
and H as an integral over the source currents in eqs. 
(UJ) makes clear that the fluctuating properties of the 
fields are directly determined by the fluctuational prop- 
erties of the source currents. In fact, if (j(r,cj)) and 
(j(r, w)j*(r',w')) are known, where the angular brackets 
symbolise the ensemble average, the corresponding aver- 
ages and correlation functions of the fields can be evalu- 
ated. Therefore, it remains to determine the stochastic 
properties of the source currents. 

From the definition of the source currents as stochastic 
Gaussian variables describing the thermal and quantum 
fluctuations inside a dielectric dissipative body, it should 
be clear that the ensemble average of the source currents 
vanishes. It follows with (j4|) that 



(E(r,w))=0 and (H(r,w))=0. 



(7) 



However, the correlations and therewith the fluctuations 
of the source currents are determined by the properties of 
the dielectric dissipative body. This connection between 
dissipation and fluctuation can be expressed in a general 
way by means of the fluctuation-dissipation theorem [l| , 
which can be applied to our problem (2) . For a homoge- 
neous, isotropic and local dielectric dissipative body this 
yields 



(j a (r,u;)j;(r>')) 

= 47tw£(cj, T)e"(u)5 a pS(r - r')S{u> - uj') 



(8) 



for the components of the source currents, where we have 
written e = e' + ie", and E(u),T) is the Einstein function 



„, mN hjj 
E(u,T) = — 



Hid 



1 



(9) 



with the usual inverse temperature /? = of the 

body. The appearance of the Einstein function is implied 
by the fluctuation-dissipation theorem, which indeed is a 
quantum mechanical operator relation, but the fluctuat- 
ing fields in the framework of fluctuating electrodynamics 
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are classical. This is a somewhat weak point in the the- 
ory, but it was recently shown 0] that a rigorous quantum 
electrodynamical treatment for dissipative media leads to 
the same correlation functions for the fields as fluctuat- 
ing electrodynamics. Hence, it seems to be justified to 
work within the framework of fluctuating electrodynam- 
ics, which should be classified as a semi-classical theory. 

Now, it is possible to evaluate the correlation functions 
of the fluctuating fields. Using eqs. (j4j together with the 
fluctuation-dissipation theorem in eq. ([8]) yields 



(E a (v,t)H (v,t)} = I dwE^T)^^-^) 



ld 3 r'\G E G Ht ) 

ct(3 



(10) 



where the f-symbol denotes the adjoint dyadic and c.c. is 
an abbreviation for the complex conjugated term. Sim- 
iliar equations can be derived for (E a (r,t)Ep(r,t)) and 
(H a (r,t)Hp(r,t)), respectively. The ensemble averages 
of the physical quantities — the Poynting vector, the 
energy density and the stress tensor — of a dielectric 
dissipative body can be evaluated with the help of these 
field correlation functions, if the material properties, i.e., 
the permittivity e(oj), and the temperature of the body 
are known. Furthermore, it is necessary to calculate the 
dyadic Green's functions for the given electrodynamical 
problem, in which the geometry of the problem is incor- 
porated. 

Radiative heat transfer between two slabs 

In principle, we are now able to calculate the radia- 
tive heat transfer between dielectric bodies of arbitrary 
shape and volume as long these bodies are homogeneous, 
isotropic, and local dissipative dielectrics, and can be de- 
scribed macroscopically. Unfortunately, the calculation 
of the dyadic Green's function for the most geometries is 
rather complicated. Therefore, we will consider only the 
radiative heat transfer between two bodies of the sim- 
plest possible geometry, i.e., two semi-infinite dielectric 
bodies separated by a vacuum gap of width d (see Fig. 

ED. 

The solution to this problem is well-known: the first 
derivation in the framework of Rytov's fluctuating elec- 
trodynamics was given by Polder and van Hove [4j . Levin 
et al. [5| solved the problem with the surface impedance 
method of Leontovich, and Loomis and Maris |6| reinves- 
tigated the problem some years ago. The radiative heat 
transfer between the two slabs can be stated as 141 



1 f°° 

(S) = -^J du [E(u,n) - E{uj,T 2 )} 

x / dk ± k ± {TP +T| 2 } 
Jo 



(11) 



where T^ 2 and T} 2 are the transmission coefficients from 
the first body at z < to the second body at z > d for 
the TM- and TE-modes, respectively. These transmis- 
sion coefficients can be expressed by means of the usual 
Fresnel coefficients Q, which gives for the propagating 
modes, i.e., for k± < ko, 
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and for the evanescent modes, i.e., for k±_ > ko, 
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4Im(rj|)Im(r|)e- 2 T d 
|1 -r^e" 2 ^ 2 : 



(12) 



(13) 



where the transmission cocfficents for the TE-modes can 
be deduced from these relations by the substitution ||— KJL. 

The mean energy transfer per unit time and unit area 
between the two semi-infinite bodies at different temper- 
atures in eq. (|lll) consists of two parts. The temperature 
enters only into the first part, which is given by the dif- 
ference of the Einstein functions evaluated at the core- 
sponding temperatures. From this difference it becomes 
clear that the vacuum term in the Einstein functions van- 
ishes; there will be no vacuum contribution to the energy 
transfer. The second part is the integral over the lateral 
wave vector k± which counts the modes contributing to 
the energy transfer. It can be shown that the above trans- 
mission coefficients are always smaller than 1. Moreover, 
it is possible to retrieve the energy transfer between two 
black bodies, which by definition have zero Fresnel coeffi- 
cients. It follows that the transmission coefficient for the 
evanescent modes (|13|) becomes zero and the transmission 
coefficient for the propagation modes (| 12[) becomes 1. In- 
serting these transmission coefficients in eq. (fTTj) yields 
the Stcfan-Bolzmann law stated in eq. ©. 

From eq. (fTTj) together with the transmission cocffi- 
cents it is possible to study the radiative heat transfer 
between the two semi-infinite bodies in detail. At the 
moment, we are only interested in the near-field limit of 
eq. (|lip . i.e., in distances kod « 1. A Taylor expan- 
sion [5j shows that in the near-field 
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and (S±) = const. 



(14) 



Therefore, the TM-mode part of the energy transfer dom- 
inates in the near-field region, S> (S±), but it is 
not clear a priori at what distance this domination be- 
gins. In fact, we will see by numerical calculations that 
the distance where the TM-mode contribution becomes 
dominant depends on the material properties of the two 
slabs. Furthermore, from the near-field limit in eq. f| 14[) 
two questions arise: i) Is a divergent radiative energy 
transfer physically reasonable? ii) Is it possible to mea- 
sure this power law? Or better, at what distance does 
the domination of the TM-mode part begin? 
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FIG. 2: Numerical results of eq. for different Drude ma- 
terials at different temperatures normalised to the black body 
value Sbb given by the Stefan-Boltzmann law in eq. ([2]). 



The first question was discussed controversially (il-[lo| . 
But it should be clear that the description based on fluc- 
tuating electrodynamics is still a macroscopic one, which 
means that at distances smaller than 10 _8 m this the- 
ory is not valid and therefore does not lead to physi- 
cally reasonable results. In fact, the source currents in 
the fluctuation-dissipation theorem in eq. ([5]) are delta- 
correlated with respect to r - r', but this can not be true 
at distances where the microscopic properties of the ma- 
terials make themselves felt. This delta-correlation of 
source currents seems to be the source of the divergent 
heat transfer at small distances, but in order to give a 
satisfactory answer to the question stated above a the- 
ory is needed that takes the finite correlations of source 
currents into account. 

With this drawback of the macroscopic theory and the 
second question in mind we will now discuss some numer- 
ical results of the heat transfer (TTTj) between two Drude 
materials For convenience we set the temperature 

of the first body to T\ = 300K and the temperature of 
the second body to T2 — OK. Moreover, we use the same 



Drude permittivity for both bodies, i.e., e% = £2- 

The numerical results for different plasma frequencies 



and d we used relatively high relaxation times and high 
plasma frequencies, so that a, b, and d are plots for good 
conductors. On the other hand, for c we used a relatively 
small relaxation time and small plasma frequency, which 
means that c is a plot for a bad conductor. 

Obviously the divergency discussed before, which is 
associated with the TM-mode part of the radiative heat 
transfer {S\\), can only be seen for c, i.e., for the bad 
conductor. There also is a divergency for good conduc- 
tors, but it appears only for distances much smaller than 
10 _9 m. Otherwise, in the region between 10 _9 m and 
10~ 7 m the radiative heat transfer for the good conductors 
a, b, and c becomes constant. As discussed before, such a 



behaviour can be attributed to the TE-mode part of the 
radiative heat transfer, (S±). As an answer to the second 
question, it follows that for bad conductors the radiative 
heat transfer in the region between 10 _9 m and 10 -7 m, 
which is accessible to measurements, is dominated by the 
TM-mode part of the radiative heat transfer, whereas for 
good conductors the TE-mode part dominates. 



SUMMARY AND OUTLOOK 

The radiative heat transfer between macroscopic di- 
electric bodies can be described by Rytov's theory of 
fluctuating electrodynamics. In the near-field the evanes- 
cent modes give the main contribution to the radiative 
heat transfer. As the numerical results indicate, for good 
conductors the TE-mode part dominates the heat trans- 
fer for experimentally accessible distances. In contrast, 
for bad conductors the TM-mode part dominates in that 
region and leads to a divergent radiative heat transfer. 
This unsatisfactory feature of the theory can be traced 
back to the delta correlation in the correlation function 
of the source currents. There is a need for a theory which 
takes a finite correlation length of the source currents into 
account, and which should lead to a finite heat transfer 
at all distances. 

Experiments now being prepared in Oldenburg should 
answer some questions and thus provide a basis for such 
a theory. Some of these questions are: At what distance 
does the macroscopic theory fail? What is the finite value 
of (S) for d — > 0? Is there an appropriate theory describ- 
ing the sensor-sample geometry? 
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